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We report on systematic numerical study of carrier multiplication (CM) processes in spheri- 
cally symmetric nanocrystal (NC) and bulk forms of PbSe and PbS representing the test bed for 
understanding basic aspects of CM dynamics. The adopted numerical method integrates our previ- 
ously developed Interband Exciton Scattering Model (IESM) and the effective mass based electronic 
structure model for the lead chalocogenide semiconductors. The analysis of the IESM predicted CM 
pathways shows complete lack of the pathway interference during the biexciton photogeneration. 
This allows us to interpret the biexciton photogeneration as a single impact ionization (II) event 
and to explain major contribution of the multiple II events during the phonon-induced population 
relaxation into the total quantum efficiency (QE). We investigate the role of quantum confinement 
on QE, and find that the reduction in the biexciton density of states (DOS) overruns weak en- 
hancement of the Coulomb interactions leading to lower QE values in NCs as compared to the bulk 
on the absolute photon energy scale. However, represented on the photon energy scale normalized 
by corresponding band gap energies the trend in QE is opposite demonstrating the advantage of 
NCs for photovoltaic applications. Comparison to experiment allows us to interpret the observed 
features and to validate the applicability range of our model. Modeling of QE as a function of pulse 
duration shows weak dependence for the gaussian pulses. Finally, comparison of the key quantities 
determining QE in PbSe and PbS demonstrates the enhancement of II rate in the latter materials. 
However, the fast phonon-induced carrier relaxation processes in PbS can lead to the experimentally 
observed reduction in QE in NCs as compared to PbSe NCs. 

PACS numbers: 



I. INTRODUCTION 

An ability to produce more then one electron-hole pairs 
(excitons) per single absorbed photon in semiconduc- 
tor materials can potentially enhance the solar energy 
conversion efficiency in photovoltaic devices beyond the 
fundamental Shockley-Queisser limitP^In the literature 
this process is referred to as Carrier Multiplication (CM) 
or multiple exciton generation, and has been extensively 
studied in the bulk semiconductors for decadesP^ Re- 
cent ultrafast optical studies of CM in semiconductor 
nanocrystals (NC) reported significantly higher QE and 
much lower activation energy threshold (AET) as com- 
pared to the bulk materials on the photon e nergy scale 
normalized by corresponding band gap valuesPSJIIH 

These reports stimulated extensive experimental and 
theoretical studies of the CM in the nanostructured ma- 
terials in which the exciton dynamics is affected by the 
quantum confinement. It was initially expected that the 
quantum confinement enhancement of the Coulomb in- 
teractions between the carriers could result in more ef- 
ficient multi-exciton production. Furthermore, a relax- 
ation of the quasi-momentum conservation constrains by 
breaking the translational symmetry should open addi- 
tional pathways for CM and reduce the AET. Finally, 
the presence of the phonon bottleneck in NCs should 
slow down the intraband phonon-assistcd cooling and fur- 
ther increase the QE. However, the subsequent reports 
claimed significan tly lower QE and even the absence of 
the CM in NCs {23HS3 This controversy possibly rises from 



experimental inaccuracy, 26 28 sample-to-sample variation 
in surface preparation) 29 " 33 ! and contrib ution from ex- 
traneous effects such as photochargingP^Ml The op- 
tical measurements of QE in bulk PbS and PbSe 
semiconductors^ show that , in fact, the bulk QE ex- 
ceeds validated values in NC a 25 1 34 1 3 ^ 2 ^ if compared on 
the absolute photon energy scaled However, due to con- 
finement induced increase in the band gap energy, E g , 
the benefits to photovoltaics is higher in NCs than in 
bulk. The controversy ca lls for development of new sen- 
sitive spectros copic toolsJ^EH alternative photo-current 
measurements J321i3] anc j reassessment of the quantum 
confinement role in the CM dynamics based on a unified 
theoretical modelP 1 ^ 1 

Accepted theoretical models of CM are based on 
the many-body Coulomb interactions which correspond 
to the valence-conduction band transitions conserving 
the total charge but not the number of electrons and 
holespSH^ In the exciton picture, this represents the 
transitions between the exciton bands of different multi- 
plicity. Hence, we will refer to these transitions below as 
the interband exciton transitions, and to the correspond- 
ing Coulomb interactions as the interband Coulomb in- 
teractions. Within adopted nomenclature, the intraband 
Coulomb interactions (i.e., restricted to an exciton band 
of a certain multiplicity) conserve the total number of 
electrons and holes and determine the multi-exciton in- 
teraction (e.g., binding) energy. 

The simplest theoretical model initially developed to 
explain CM in bulk materials is the impact ionization (II) 
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modelf^ It treats CM dynamics following high energy 
photoexcitation as the lowest order interband transitions 
whose rate is given by the Fermi's Golden Rule. This 
rate depends on the interband Coulomb matrix elements 
and the final multi ple-carriers (e.g., biexciton) density of 
states (DOS)p21S2l xhe model has been further applied 
to interpret the CM dynamics in NCs along with various 
methods for the electronic structure calculations. 

Franceschetti, An and Zunger performed atomistic 
pseudopotential calculations of the II and Auger recom- 
bination rates in a PbSe NC. Assuming the absence of 
the quasi-momentum constraint and rapid growth of the 
biexciton DOS, they estimated low value (~ 2.2E g ) of 
AET defining this quantity as the energy above which the 
II rate exceeds the Auger recombination oneP Rabani 
and Baer performed screened semiempirical pseudopo- 
tential calculations of the II rates in CdSe and InAs and 
found significant reduction of the QE with the NCs size 
increase. Observed trend makes CM already inefficient 
in NCs of diameter ~ 3 nm. The observation has been 
rationalized by strong size dependence of the interband 
Coulomb interactions and the trion DOS behaviorPS 

Using the semiempirical tight binding model Allan and 
Delerue performed exte nsive calculations for PbSe and 
PbSe, InAs and Si NCsPH They found that although 
the Coulomb interactions are enhanced by the quantum 
confinement, the quantum-confinement-induced reduc- 
tion in the biexciton DOS facilitates the reduction of the 
II rate and subsequently of QE. The comparison of the II 
rates calculated in PbSe bulk and PbSe NCs shows that 
the II rates in the NCs does not exceed the bulk oneP 
Based on these calculations, the authors of Ref . [315 argue 
that the experimentally observed drop of QE in NCs can 
be attributed to the dominant effect of the reduced car- 
rier's DOS. Lack of the phonon bottleneck, leads to the 
rapid intraband phonon-assisted relaxation which further 
reduces the QEPHH 

Schaller, Agranovich, and Klimov, first pointed out 
that in addition to II, the direct photogencration of bicx- 
citons can take place in NCs through "virtual-exciton 
channel".^ Subsequently, Rupasov and Klimov sug- 
gested that additional contribution to photogeneration 
QE can rise from the "vitual biexciton channel" Us- 
ing the effective mass model for the electronic structure 
of PbS and PbSe NCs due to Kang and WiseJ^ the au- 
thors estimated the photogeneration QE and argued that 
this process dominates the CM dynamics in NCs. Re- 
cently, Silvestri and Agranovich using the same model, 
performed detailed calculations for relatively small radii 
(specifically for 1.95 nm and 3 nm) PbSe NCs P They 
concluded that the contribution of the photogeneration 
processes in NCs is much weaker than it was claimed be- 
fore, and clarified that the overestimated QE results from 
disregarding the effect of selection rules for the interband 
Coulomb matrix elements and the oscillator strength fac- 
tors weighting optically allowed transitions. Here, we 
demonstrate that the absence of the interference between 
the photogeneration pathways and the presence of small 



size dispersion further reduces the photogeneration QE. 

Another model describing the coherent photogenera- 
tion of biexcitons from the resonant exciton states ini- 
tially proposed by Shabaev, Efros, and Nozik considered 
idealized five level systemp^ and was subsequently re- 
fined to account for the biexciton DOS effect Using the 
effective mass k • p Hamiltonian, the authors of Ref. [59"] 
performed calculations for small radii (i.e., 2 nm and 
3 nm) PbSe NCs, and found that dense biexciton DOS 
leads to the vanishing coherent oscillations. They also 
pointed out that efficient CM can be expected in consid- 
ered small NCs. The drawback of the coherent superpo- 
sition model is the lack of the pure-dephasing effects and 
the inhomogeneous broadening which as we demonstrate 
here play important role in the CM dynamics reducing 
the QE. 

Finally, the ah initio calculations on small clusters 
(< 1 nm) demonstrated the role of strong Coulomb corre- 
lations and fast excit on-bi exciton dephasing in multiexci- 
ton photogeneration} 60 ! 61 ! and suggested contributions of 
the phonon-assisted Auger processes to CMp^The atom- 
istic calculations mostly focused on II dynamics occurring 
during the population relaxation are limited to the small 
diameter (< 3 nm) NCs. However, reported experimen- 
tal studies consider larger NCs with the diameter varying 
up to 10 nm. To close this gap and to perform compari- 
son with bulk, an extrapolation procedure combined with 
atomistic calculations has been usedl^The photogenera- 
tion dynamics has been investigated using effective mass 
models for larger but still small diameter (^4 — 6 nm) 
NCs. 

Recently, we have published a letter summarizing the 
results of our numerical investigation on the direct pho- 
togeneration and population relaxation processes con- 
tributing to CM in a broad diameter range PbSe NCs 
and in PbSe bulkP The calculations employ our earlier 
proposed Interband Exciton Scattering Model (IESM)P 
which recovers the models discussed above as limiting 
cases and which we further parameterized using the ef- 
fective mass, k • p, electronic structure model proposed 
by Kang and WiseP^^ In the letter, we argue that in 
both cases of the photogencration and population relax- 
ation the II is the main mechanisms of CM. This explains 
weak contribution of the direct photogeneration to the 
total quantum QE. Analyzing the scaling of the total 
QE with the NC size, we found that QE in NCs plot- 
ted on absolute energy scale does not exceed that in the 
bulk. This is in agreement with reported experimental 
data and some theoretical studies, and confirms that the 
quantum-confinement-induced reduction in the biexciton 
DOS makes a dominant contribution to QE. 

Present paper discusses in great details the theoretical 
method and provides extended analysis of the numerical 
results summarized in Ref. [64 j. In addition, we con- 
sider the effect of the optical pump pulse duration on 
QE, and provide a comparative analysis of CM in PbS 
and PbSe materials. The paper is organized as follows. 
In Sec. [Hj we review the weak Coulomb interaction limit 
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of the IESM used for numerical calculations, introduce 
convenient quasicontinuous representation, and further 
discuss the details of the model parameterization. In 
Sec. |III| an extensive analysis of our numerical results on 
PbSe NCs is given. Comparison of the key CM features 
in PbSe and PbS is performed in Sec. |IV| Sec. [V] con- 
cludes the paper by a discussion on the limitations of the 
adopted model, connection with experiment, and possi- 
ble improvements of QE in NCs. 



II. THEORETICAL MODEL AND NUMERICAL 
PARAMETERIZATION 

The central quantity describing CM efficiency is QE 
which can be calculated as 



QE 



2N X 



N x 



N x 



N x 



(1) 



Here N x and N xx are the total exciton and biexciton 
populations, respectively, produced by single absorbed 
photon in the limit of vanishing pump fluencies. They 
depend on the delay time measured from the center of 
the pump pulse and allow one to determine both the QE 
due to the photogeneration event and the total QE after 
the population relaxation. A complimentary quantity, 
often used to describe CM yield, is the biexciton quantum 
yield, r\ xx = QE — 1, defined as a number of biexcitons 
produced per single absorbed photon. 

In this section, we describe the weak Coulomb limit 
of the IESM 47 which we use to calculate the time evo- 
lution of the exciton and biexciton populations entering 
Eq. ([I]) . We also introduce the quasicontinuous frequency 
representation which is convenient to analyze the quan- 
tum confinement signatures and the size scaling of the 
quantities determining the QE in transition from NC to 
the bulk limit. Finally, we discuss the model parameter- 
ization for the numerical calculations of QE in PbSe and 
PbS materials. 



is the pulse self-convolution function depending on the 
ground state to exciton state transition frequency, u> x , 
and dephasing rate, j x . The pump pulse is characterized 
by the envelope amplitude, £ pm (t), giving the temporal 
profile, and the central frequency, Lu pm - In our calcula- 
tions, we use the Gaussian form of the pulse envelope 
function, £ pm {t) = £ P m,e~ t / 2t p»", with the mean am- 
plitude, £pm, and pulse duration, r pm . For the pulses 
longer than typical dephasing time, r pm 3> I/7, the 
continuous wave (cw) limit is recovered and the pulse 
self-convolution function becomes proportional to the 
Lorentzian line shape function. 

To be consistent with the calculations of the photo- 
generated biexciton population appearing in the second 
order interband Coulomb coupling, we account for the 
first and second order corrections to the exciton popu- 
lation and disregard the effect of the exciton coherences 
both introduced in Ref . [37] . The negligible contribution 
of the latter quantities has been validated by our direct 
numerical calculations. 

The photogencrated population of the k-th. biexciton 
state is33 
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where fi^f is the intraband transition dipole between the 
biexciton states \xxk) and \xx{), and 
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A. Biexciton Photogeneration and Population 
Relaxation Model 

In the weak Coulomb interaction limit, the leading con- 
tribution to the population of the a-th exciton state gen- 
erated by the pump pulse is given by 

< i0) = | M *o| 2 ^a-^ Pm ), (2) 

where [i x is the transition dipole moment between the 
ground state, \x ), (i.e., filled valence band) and an exci- 
ton state \x a ), and 

c\ /"OO /"OO 

I(u% - Upm ) = ¥ J dt J dtie-<^ (3) 
x cos[(w£ -Upm)h]£pm(t )£ pm {t 



describes the transition amplitude associated with the 
single (Born) scattering event between exciton sate, \x a ), 
and biexciton state, \xxk)- The quantity A^q' x , also en- 
tering Eq. Q, is the transition amplitude between the 
ground state, \xq) (ui x = 0), and biexciton state, \xxi). 
V^ x a ' x is the interband Coulomb matrix element in the 
exciton representation giving rise to the interband scat- 
tering, and is the interband dephasing rate. The 
pulse self-convolution function, T-iuf^ — w pm ), entering 
Eq. Q contains u)f x and ^ xx instead of oj^ and j x , re- 
spectively. According to our calculations the coherence 
contributions to the biexciton population^! are also neg- 
ligible. 

Our numerical calculations show that the rapid sign- 
variation of the interfering terms in Eq. Q leads to their 
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FIG. 1: Non-interfering biexciton photogeneration pathways associated with the weak Coulomb limit, (a) Indirect biexciton 
photogeneration via exciton states, (b) Direct biexciton photogeneration via virtual exciton states, (c) Direct biexciton 
photogeneration via biexciton states. (See details in the text.) 



cancelation allowing us to write Eq. Q as 

I 2 
xx(2) _ \aXX,x x J(,.x_.. \ 

2 
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The resulting three non-interfering photogeneration 
pathways are illustrated in Fig. [T] 

The first and the second pathways, given by the first 
and the second terms in Eq. ^ , are shown in panels (a) 
and (b), respectively. These terms describe redistribu- 
tion of the exciton oscillator strengths (~ I Mora 1 2 ) between 
exciton and biexciton bands mediated by the interband 
Born scattering (A x,xx ). First pathway (panel (a)), in- 
volves the ground-state-to-exciton resonant optical tran- 
sition and further scattering to final biexciton states dis- 
tributed around hcj pm according to the non-zero compo- 
nents of A x,xx . We will refer to this process as the indi- 
rect biexciton photogeneration throughout this paper.^ 
Panel (b) illustrates the second pathway, where the exci- 
ton is virtual and final biexciton state is in resonance with 
the optical pulse. Thus, this pathway will be refereed to 
as the direct biexciton photogeneration via virtual exciton 
states!^ The third pathway (panel (c)) consists of Born 
scattering between the ground state and biexciton states 
(A 0,xx ) stabilized by the intraband dipole transition, fi xx . 
Accordingly, the final biexciton state is in resonance with 
the optical pulse. This process, refereed to as the direct 
biexciton photogeneration via biexciton states]^ becomes 
prohibited in the bulk limit by the momentum conserva- 
tion constraint for optical valence-conduction band tran- 
sitions. 

Equation ^ with the higher order corrections and 
Eqs. Q-(|6]) fully describe the exciton and biexciton pop- 
ulations prepared by the pump pulse. We use them for 



the numerical evaluation of the photogenerated exciton, 
N x = J2a n a> and biexciton, N xx = J2k n T > popula- 
tions and QE (Eq. ([!])). They also provide the initial 
conditions to model the population relaxation dynamics 
using a set of rate equations Using this computational 
approach, our goal is to clarify the effect of quantum con- 
finement on QE in transition from NCs to the bulk limit. 
Since the CM dynamics in NCs takes place in the energy 
region characterized by high electron and hole DOS, we 
expect that the CM dynamics in NCs and in the bulk 
should have common features. Thus, we intend to see 
how strongly these features are affected by the presence 
of the confinement potential. 

Quantitatively, we are going to look at the interplay 
between the size scalings of the interband Coulomb inter- 
action and exciton/biexciton DOS determining the QE 
variation in transition from NC to the bulk limit. First, 
we define the bulk limit as the thermodynamic limit: 
V — > oo, V/v — > co, and v = const., where V is a crystal 
volume, v is the unit cell volume, and the ratio V/v gives 
the number of unit cells in the crystal. Then Eqs. ^ 
and j6j|, and the quantities determining the population 
relaxation should be represented in such a form that the 
effect of the interband Coulomb interactions and exci- 
ton/biexciton DOS are clearly distinguished. This can 
be achieved by using the quasicontinuous energy repre- 
sentation. Since some quantities of interest in the bulk 
limit have volume scaling, it is convenient to introduce as- 
sociated intensive (i.e., volume- independent in the bulk 
limit) variables. Their deviation from the well defined 
bulk values will provide us with the convenient measure 
of the quantum confinement effects. 

We start with the exciton and biexciton populations 
which in the quasicontinuous energy representation read 
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respectively. Both of them have linear scaling with the 
volume in the bulk limit. Therefore, we eliminate the 
latter scaling by multiplying the total populations with 
the dimensionless prefactors (v/V) and end up with the 
following intensive quantities 



N x = 



V 



n x (uj)duj, 



— I n xx (u))du. 
v Jo 



(9) 
(10) 



Although the prefactor v/V cancels out in the expression 
for QE (Eq. ([!])), indicating that the latter quantity is 
indeed intensive, it is convenient to keep it in Eqs. ^ 
and ( |To| for consistency. 

According to Appendix [Aj we define the intensive ex- 
citon and biexciton DOS as 



V 



(^) 4 x> 



I , xx \ 
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respectively. The associated optically allowed exciton and 
joint biexciton DOS read 



px(u) 

P~xx(ui,U) 2 ) 
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respectively. They carry information on the optical se- 
lection rules reducing the number of states participating 
in the photogeneration process. 

Our central quantity is the effective Coulomb term de- 
fined as the r.m.s. of the interband Coulomb matrix ele- 
ments which couple states within the frequency intervals 
[u>l,w\ + du>i] and [ll^,^ + du>2], 



V~x,xx / \ 
e j f (0Ji,UJ2) 



El rrx,xx I 
| *a,m | 



(15) 
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Related effective Coulomb term, coupling the ground 
state and biexciton states, can also be defined as 



T rxx 
V eff 
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(16) 



As we discuss in App end ix [A"} the volume prefactor 
(V/v) 2 [(V/v) 3 / 2 } in Eq. (JTH]) [Eq. |l6|] corresponds to a 
finite effective Coulomb value in the bulk limit. There- 
fore, the size scaling of such defined interaction with the 
NC diameter, d, provides quantitative measure of the 



quantum confinement effects. In general, a deviation 
from the bulk value for the effective Coulomb reflects 
the net result of the scaling of the Coulomb matrix ele- 
ments with d, relaxation of the momentum conservation 
constraints, and the appearance of new selection rules 
associated with symmetry of the confinement potential. 

Assuming the cw excitation, the exciton population 
produced by the pump pulse is simply proportional to 



the corresponding optically allowed DOS (Eq. ( 13 ) ) , 

N x (u> pm ) = Ap x (u) pm ), 



(17) 



where A = 2n 3 / 2 T pm £$ l /ti 2 . The biexciton population 
as a function of the pump frequency in the adopted rep- 
resentation reads 



Nxx (jdpm ) 
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As desired, Eq. ( 18 ) is volume independent and clearly 



distinguish the contributions from the effective Coulomb 
interactions and variously defined DOS. Accordingly, 
Eqs. (17) and (18) provide central expressions for the 



interpretation of the numerical simulation results on the 
photogeneration QE as discussed in the subsequent sec- 
tion. 

The population relaxation dynamics is described by 
the following set of kinetic equations 

h x (uj) = -k n (uj)n x (uj) + k AR (oj)n xx (uj) (19) 

duj' T x (uf ' ,u>)n x (u) 



du>' T x {u,(J)n x {tJ), 
h xx (u) = k n (u)n x (ui) - k A R{u)n xx (uj) (20) 
dw' T xx (lj' ,u)n xx (u)) 

dJ T xx (ui,uj')n xx (uj'). 

These equations contain both the II and Auger recombi- 
nation rates 



kAR{u) 



2tt 

-^[V x f f x (u)] 2 p xx (u), 



(21) 



^ (^) Ve/f MfM"), (22) 



h 2 



respectively. Here the shorthand notation V x j xx (ui) 
stands for the diagonal component of the effective 
Coulomb term, i.e., V x j xx (iS) = V x J xx (oj,ui). According 
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to Eq. (21 1, the II rate is an intensive variable. In con- 



trast, the Auger recombination rate (Eq. (22)) vanishes 
in the the bulk limit as d~ 6 . 

In the case where the II processes dominate the CM dy- 
namics, the ratio of the II to Auger recombination rates 



kniuj) 



* Alt 



(23) 



should determine QE. Since this ratio is proportional to 
the ratio of the corresponding DOS, it has been proposed 
as a selection criterion for the materials showing efficient 
CM.HH However, Eq. ( 23 ) clearly shows that besides the 



material-specific signatures given by the ratio of the in- 
tensive DOS, it also contains the volume scaling factor. 
Therefore, we argue that Eq. ( 23 ) can only be used for a 



material selection criterion, after the volume prefactor is 
eliminated. 

The intraband relaxation rates 



in Eqs. (19) 



and (20) describe the phonon- assisted cooling. In the ab- 



sence of the phonon bottleneck and in the region of high 
exciton DOS, it is expected that single-phonon processes 
dominate the exciton and biexciton intraband relaxation. 
Accordingly, we ca lculat e these quantities by using the 
following expressio n 1 66 * 67 ! 



Ti(uj\ uj) = sign(u/ — lu) 



2J i (|w'-w|) 



exp[H(u>' — uj)/kBT] — 1 ' 



(24) 



where i ksT is thermal energy, and Ji(Aw) is 

the phonon spectral density approximated by the ohmic 
form with exponential cutoff 
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Here, the adjustable parameters are electron-phonon cou- 
pling, Ai, with the constraint X xx = 2X X and the phonon 
frequency cutoff, oj c . 



Eqs. (19)-(25) represent the quasicontinuous represen- 



tation of the discrete kinetic equations introduced in 
Ref. 07] ■ We found that the former representation is 
more suitable for numerical integration. The solution of 
these equations with the initial conditions determined by 
the photogenerated populations (Eqs. ([2])-(|6])) provide 
closed computational scheme that we use to determine 
the total QE. Next, we discuss the parameterization of 
the introduced model and some details on the adopted 
numerical techniques. 



B. Model Parameterization 

An accurate knowledge of single electron and hole wave 
functions is required to construct the exciton and biex- 
citon states and to evaluate transition dipoles and the 
interband Coulomb matrix elements. For this purpose, 
we adopt the effective mass k • p formalism, originally de- 
veloped by Mitchell and WalliiP^ and DimmoclpS for the 



bulk lead chalocogenide semiconductors, and further ad- 
vanced by Kang and Wise^ for the spherically symmet- 
ric NCs. The formalism is based on a four-band envelope 
function model explicitly taking into account spin-orbit 
interaction between valence and conduction bands. In 
considered PbSe and PbS materials, the band structure 
anisotropy is small, and therefore, neglected in our cal- 
culations. 

An electron and hole wave function obtained within 
this formalism reads 



Mr) = £ F; n (r)u m {v). 



(26) 



m— 1 



Here u m (r) is m-th component of the bulk Bloch wave 
function associated with the band-edge states i n L- v alley 
whose index, m — 1, . . . , 4, denotes the bandsP^! The 
envelope eigenfunctions i^(r) and the corresponding 
eigenenergies, Hu>i, are found by solving the k • p Hamil- 
tonian eigenvalue problem with the infinite wall bound- 
ary condition, Fi(\r\ = R) = 0, at the surface of a spheri- 
cal NC of radius Rp^If frwi > (hui < 0) we identify the 
state as a conduction band electron (valance band hole) 
state. The eigenstate index represents a set of quantum 
numbers, i = {n, n,j, m}, such as primary quantum num- 
ber describing number of the wave function nods, parity, 
total angular momentum, and its projection, respectively. 
In what follows, we will refer to the single particle states, 
described by Eq. (26), as the Kang- Wise (KW) states. 



In the bulk limit, we set the envelope functions, F^ir), 
to the plane waves. 

A natural way to introduce the exciton and biexci- 
ton states is to use the second quantization representa- 
tion within the basis of the KW-states (Eq. (Bl )). Tak- 
ing into account that the electron-hole Coulomb inter- 
actions are weak compared to their kinetic energies (i.e., 
strong confinement regime) , we introduce exciton and bi- 
exciton states as the following configurations of uncorre- 
cted electron-hole pairs 



|x ) = c\dl\x ), 



\xx k ) 



d\dl\x ). 



(27) 
(28) 



As introduced above, the ground state, \xq), is attributed 
to the fully filled valence band. An exciton (biexciton) 
energy in terms of the electron, fkOp, and the hole fuv^ 
energies is simply: Hlu x = hcu^ — ftu^ (fiw% x — Huj^ + 
huj e q — huj^ — huj^ ) . Note that the exciton index a = {p, q} 
(biexciton index k = {pq, rs}) is a collection of two (four) 
sets of KW quantum numbers. 

The expressions used to calculate the dipole moments 
/U^q and /Ufef, and the interband Coulomb matrix ele- 
ments, V x < xx , in Eqs. @, {5}, ([l5|, and (16} using 
associated matrix elements in the KW basis are given 
in Appendix [5] Numerical calculations of the interband 
Coulomb KW matrix elements are performed by using 
the multipole expansion in terms of the Clebsch-Gordan 
coefficients. The dephasing rates for the high excited 
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states are not available from experiment, and in our cal- 
culations we use some average value of ftf = fvy xx /2 = 
tvy x ' xx /2 = 50 meV. This value is con sistent with the 
calculations performed on PbSe clusters!- 60 * 70 * 71 * 

The numerical calculations of the photogenerated pop- 
ulations have been performed according to Eqs. ([2])-j5]) 
including the first and second order corrections to Eq. pF. 
To account for the degeneracy of the exciton and biexci- 
ton states associated with the four equivalent L-valleys, 
we multiply the r.h.s. of Eq. ([2| and Eq. Q by factors 4 
and 16, respectively. This assumes that the inter-valley 
Coulomb scattering processes are negligible. In most cal- 
culations, the pump pulse duration is set to r pm = 50 fs, 
which is typical value used in the experimental studies.^ 
The number of exciton and biexciton states to handle 
computationally is extremely high. Even for PbSe NCs 
of moderate size (<i ~ 5 nm), the number of the biex- 
citon states with energies below 4E g is ~ 10 5 . In the 
bulk limit, the number of the exciton and and biexciton 
states becomes infinite rendering the direct summation 
over these states in Eq. Q impossible. Therefore, in 
both cases of finite-size NCs and the bulk, we use the 
Monte Carlo sampling which allows us to consider sig- 
nificantly larger NCs than it was possible in the work of 
Silvestri and AgranovichP^ 

To perform numerical calculations of the popula- 
tion relaxation processes, we first evaluated the in- 
terband Auger recombination and II rates given by 
Eqs. (21|-(22|. For this purpose, we smeared the delta- 

and 



functions in the effective Coulomb term (Eq. ( jl? 
exciton/biexciton DOS (Eqs. (Ill and Jl2|) with the 
Lorentzian profiles containing the corresponding dephas- 
ing rates and used the Monte Carlo procedure to sum over 
the exciton and biexciton states. To evaluate the phonon- 
induced intraband relaxation rates (Eqs. (24) and (25)), 
we set the room temperature value for fcgT = 25 meV, 
the cutoff energy to hu> c = 50 meV. The electron-phonon 
coupling X x — X xx /2 was fit to reproduce the experimen- 
tally observed intraband relaxation time in the range of 
0.5 < Tp h < 5 pspfiESl The kinetic equations (Eqs. (l9| 



and (20)) have been further numerically integrated on 
the energy grid. All calculated observables are averaged 
over ensemble of NCs with the Gaussian size distribution 
with the standard deviation a = 5%. 



III. INTERPLAY OF CM PATHWAYS IN PbSe 

In this section, we discuss the results of our numer- 
ical calculations of the photogencration and population 
relaxation processes in NC and bulk PbSe to clarify the 
effect of the CM pathways interplay and quantum con- 
finement on QE. As we already mentioned, the QE is 
determined by the interplay between the DOS and the 
effective Coulomb scalings. Therefore, to set the stage 
for our analysis we first consider these quantities sepa- 
rately. 

Calculated energy dependence of the exciton and biex- 
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FIG. 2: Energy dependence of (a) exciton and (b) biexciton 
DOS calculated for various diameter, d, PbSe NCs and the 
bulk, (c) The diagonal component of the effective Coulomb 
term in PbSe NCs and the bulk. The inset shows actual 
interband Coulomb interaction in NCs defined as the effective 



Coulomb term with the volume prefactor added, i.e.. 

(v/v) 2 v e y f x . 



v, 



c 



citon DOS (Eqs. ( fTTj ) and (12)) for different diameter 
PbSe NCs and the bulk limit are shown in panels (a) 
and (b) of Fig. |2j respectively. Obviously, the quantum 
confinement leads to the DOS reduction in NCs com- 
pared to the bulk. The energy dependence of the DOS 
follows power law. Specifically for the bulk limit, we 
find that p x (oj) ~ (w — w g ) 2 ' 2 and p xx (w) ~ {oj — 2w g ) 6 
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FIG. 3: (a) The pump energy dependence of the photogener- 
ation QE for various diameter, d, NC and bulk PbSe plotted 
on the absolute pump energy scale. Solid lines represent the 
results obtained from numerically exact (Eqs. |2|-([5}) cal- 
culations, and the dash indicates approximate calculations 
(Eq. (30l) for d=5,8,15 nm NCs and the bulk using the effec- 
tive dephasing rates tv~f & ff = 50,40,55, and 45 meV, respec- 
tively, (b) Solid lines from (a) plotted on the unitless pump 
energy scale (i.e., normalized by corresponding NCs and bulk 
band gap energies, E g ). 



with hujg = 0.28 eV reflecting contributions of the non- 
parabolic regions of the electron/hole band structures. 
Even in the high DOS regions where the bulk behavior 
is expected in NCs, the DOS converges rather slowly to 
the limiting bulk values as the NC diameter increases. 

The effective Coulomb as a function of energy is plotted 
in Fig. |2j (c) and follows power law behavior, V^ x (uj) ~ 
w~ 5 . According to the plot, the effective Coulomb in- 
teraction is only weakly enhanced by quantum confine- 
ment, except for small, d — 5 nm, NCs. Specifically, the 
confinement results in a factor two enhancement in NCs 
with d = 8 nm as compared to the bulk. Therefore, the 
scaling of the actual interband Coulomb interaction, i.e. 
V c = {v/V) 2 V^f£, in NCs (inset to Fig. g (c)) in the 
region of high DOS should be dominated by the vol 
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pm 



FIG. 4: Contributions of the photogeneration pathways to the 
biexciton quantum yield calculated for d — 15 nm PbSe NC. 
Black and red lines show the contributions of the two com- 
ponents associated with the first two pathways schematically 
shown in panels (a) and (b) of Fig. [I] respectively. The green 
line shows the contribution of the third pathway illustrated in 
panel (c) of Fig. [I] The inset presents the relative contribution 
of the third pathway to the total photogeneration quantum 
yield in various size NCs. 



prefactor, i.e., ~ d~ 6 . However, we remind that the vol- 
ume prefactor responsible for such dramatic scaling does 
not enter the quantities determining the QE except the 
Auger recombination rate (Eq. (17)-(22)). 



A. Photogeneration QE 

The photogeneration QE as a function of the abso- 
lute photon energy in NC and bulk PbSe is shown in 
Fig. [3^a). According to the plot, the QE at fixed pho- 
ton energy monotonically increases with the NC diameter 
but does not exceed the bulk value. In contrast, the QE 
plotted on the unitless photon energy scale (Fig. |3jb)) 
shows the opposite trend useful for photovoltaic applica- 
tions. However, the absolute energy scale representation 
is more suitable to understand the physical mechanisms 
of photogeneration and, therefor, we continue our anal- 
ysis using this scale. First, we compare our calculations 
for small, d = 5 nm, NC with the calculations reported 
by Silvestri and Agranovich 57 for d — 6 nm PbSe NC. 
In Ref. [57], the size dispersion of the NCs is set to 2% 
resulting in the sharp peaks reaching QE w 1.2. The in- 
crease in the dispersion to realistic value of 5% (red solid 
line in Fig. [3^a)) washes out the peaks and reduces the 
QE by one order of magnitude, i.e. to QE ~ 1.01. 

The contributions of different pathways (Fig. [I]) to the 
biexciton photogeneration quantum yield, rj xx = QE — 1, 
are depicted in Fig. [4] First of all, we point out that for 
b" the volume the energy region fiw pm > 1.5 eV in which the biexciton 
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tion via biexciton states (green line) has negligibly small 
contribution compared to the other terms. The inset 
shows that the relative contribution (i.e., normalized per 
total photogeneration QE) of the latter pathway for all 
considered NC sizes is small and has tendency to de- 
crease with the NC size growth and to vanish in the bulk 
limit. This contribution becomes small due to the strong 
off-resonant nature of the ground state to biexciton state 



effective Coulomb term (Eq. ( 16 1) entering the third term 
in Eq. (fill. 



The small contribution of the last pathway leads us 
to an important conclusion that the photogeneration QE 
is fully determined by a single II event occurring on the 
dephasing timescale. According to Eqs. (fij) and (29) this 
QE can be approximated by the following simple expres- 



FIG. 5: Calculated photogeneration QE as a function of pump 
energy for PbSe NC of d = 15 nm and PbSe bulk. For compar- 
ison, we also show experimentally measured total QE in PbSe 
NCs^ with the diameter varied in the range of 5 < d < 8 nm 
and in the bullP^. The inset shows calculated QE for PbSe 
NC of d = 5 nm on extended pump energy scale compared to 
the experimentally measured total QE for the same size NC. 



DOS shows steep growth, the indirect biexciton photo- 
generation (black line) and direct biexciton photogener- 
ation via exciton state (red line), provide identical con- 
tributions. This behavior can be interpreted by looking 
at the related first and second terms in Eq. ( 18 1 in which 



the resonant nature of the denominator, (u — w p 
corresponds to the leading contribution of the diagonal 
component of the effective Coulomb term V^ x (uj pm ) = 
Vgjf X (oj pm , ujp m )- Hence, the integral convolutions can 



be evaluated resulting in the identical two terms whose 
net contribution to the biexciton population is 



N xx {u pm ) = kll(uJpm)N x (ujp m )/~f e ff. 



(29) 



Here, ku is the II rate given by Eq. (21 1, and 7 e // is ef- 



fective interband dephasing rate. Derived expression has 
a very clear physical interpretation: The photogenerated 
biexciton population associated with the first and second 
pathways is a result of a single II event taking place on the 
dephasing timescale, l^ff a-nd following optical prepara- 
tion of the exciton states^l 

The direct biexciton photogeneration via exciton states 
is manifested by the direct dependence of the biexciton 
DOS on uj r , 



u pm in the second term of Eq. ( 18 ). As a result, 
for hoj pm < 1.5 eV, this contribution (redfine in Fig. Q) 

1.5 eV can be iden- 



pra 



drops quickly. Accordingly, hui 
tihed as the AET for the photogeneration processes in 
d — 15 nm NCs. In contrast, the indirect biexciton pho- 
togeneration (black line in Fig. Q) is resonant at exciton 
states, and decreases relatively slow for Hu) pm < 1.5 eV. 
According to the first term in Eq. ( |18[ ), this behavior re- 
flects the behavior of the off-diagonal effective Coulomb 
component V^fl (u) ,ui pm ) at Huj < 1.5 eV. 

According to Fig. [4j the direct biexciton photogenera- 



QE(uj pm ) 



2kn(uj p 



"feff 



kll{u)pm) + Jeff 



KII 
leff'' 



(30) 



where ku is the II rate given by Eq. (21), and Jeff is 
effective interband dephasing rate. Notethat the optical 
selection rules do not enter Eq. (30). The approximate 



form of QE given in Eq. (30) is quite general, since the 



approximation is based on a general fact of weak inter- 
band Coulomb interaction. To verify this relation, we 
calculated the QE using Eqs. pi] ) and (30), and com- 
pare the results (dashed lines) in Fig. [3jayT^ In general, 
the dashed lines well reproduce the trends in the behav- 
ior of the associated solid lines confirming the validity 
of the approximation. Observed small discrepancies are 
due to the numerical noise, phenomcnological origin of 
the effective dephasing rate neglecting its frequency de- 
pendence, and neglect of the correlations between the 
effective Coulomb and DOS fluctuations during the NC 
ensemble averaging. 

Both direct biexciton photogeneration pathways 
via e xciton a nd biexciton states have been studied 
beforeP ^ * 55 * 57 * The indirect biexciton generation as we 
demonstrated above is also important and should not 
be omit ted. Initial studies of the photogeneration 
pathway d 15 * 55 ! also based on the KW parameterization 
reported significantly larger QE compared to our results 
and the results report in Ref . [57] . The cause of the over- 
estimate is the disregard of the Coulomb coupling se- 
lection rules and the oscillator strengths factors weight- 
ing optically allowed transitions.^ Next, we show that 
the interference related power-scaling of the DOS not ac- 
counted for in the previous studies further leads to sig- 
nificant reduction in the photogenerated QE. 

Provided, the terms in Eq. Q associated with the 
direct biexciton photogeneration processes via exciton 
and biexciton states are interfering constructively, one 
can show that associated biexciton populations scale as 
[p x ] 2 p X x and [p xx ] 2 Pxx, respectively. This is exactly the 



case considered in Refs. |15l55j Here, p x (p xx ) is DOS for 
the intermediate exciton (biexciton) states, and p xx the 
biexciton DOS at the pump frequency. Since, we found 
that there is no interference between the photogenera- 
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FIG. 6: Energy dependence of the calculated II rates in PbSe 
NCs of different diameter and in PbSe bulk. The inset shows 
associate d A uger recombination rates calculated for NCs us- 
ing (Eq. J23b). 



tion pathways, the first two terms in Eq. (18 1 scale as 
PxPxx-, and the last one as p xx ~ p xx p xx . The latter 
linear dependence on the intermediate exciton and biex- 
citon DOS versus former quadratic on e! 15 * 55 ! significantly 
decreases the biexciton photogeneration quantum yield 
and, in fact, allows us to interpret the photogeneration 
dynamics as a single II event (Eq. (29))! 

In Fig . IS"! w e compare the experimentally measured 
total QE 34 I 36 I with the calculated photogeneration QE. 
The comparison with experimental data for NCs is pos- 
sible only at the photon energy fixed at ~ 3.1 eV, where 
the vertical aligned diamonds represent the QE measured 
for the NCs with diameters varying in the range between 
5 nm and 8 nm. The comparison shows that for the bulk 
and for the NCs of comparable diameter the experimen- 
tally measured total QE significantly exceeds the calcu- 
lated photogeneration QE. Accordingly, we conclude that 
the biexciton photoexcitation along does not explain the 
experimentally observed QE values. Next, we investigate 
the contribution of II events to the total QE taking place 
during the phonon-assisted cooling. 



B. Effect of population relaxation and pump pulse 
duration on total QE 

Central quantity defining QE during the population 
relaxation (and as concluded above during the photogen- 
eration event as well) is the II rate (Eq. (21 )). Fig. [6] 
shows calculated energy dependence of the II rate for 
NC and bulk PbSe. According to the plot, the rate 
scales linearly away from AET. The scaling directly fol- 
lows from Eq. (21) in which one can insert the effec- 

, ,-5 



tive Coulomb, 

Pxx{u) ~ (CJ - 
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2u) g ) e , whose scalings were determined 



FIG. 7: Pump energy dependence of the total QE in PbSe. 
(a) Calculations performed for bulk PbSe using various intra- 
band relaxation times, T p h- Experimental data from Ref. |36) 
are shown for comparison, (b) Calculations for both NCs and 
bulk using r v h = 1 ps (solid lines) and t v h — 2 ps (dash). For 
comparison the experimental data from Refs. [34] and [36] are 
plotted in diamonds and dots, respectively. 



at the beginning of this section (Fig. [2]). In the lit- 
erature, quadratic and even higher power behaviors of 
the II rate are obtained near the AET. As proposed in 
Rcfs. 50 72 73 such scalings solely follow from the DOS 
behavior. However, these theories do not account for the 
long-range interband Coulomb corrections at the energies 
higher than AET. As we show in Appendix |B| the first 
non-vanishing k • p contribution to the matrix elements 
of the interband Coulomb term yields additional factor 
lu^ 1 / 2 . As a result, one gets ku(oS) ~ (to — loq) 2 /lo ~ uj 
at uj 3> wo j which is in agreement with our numerical 
calculations. 

According to Fig. [6j the II rate increases as the NC 
diameter increases but does not exceed the bulk values. 
This observation is a result of the weak effective Coulomb 
enhancement in NCs observed in Fig . [2] (c) . In agreement 
with previously reported studies,^ we find that this en- 
hancement is fully suppressed by the reduction in the 
biexciton DOS (Fig. [2 (b)). As we already pointed out, 
strong size scaling of the actual interband Coulomb in- 
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FIG. 8: Comparison of the calculated total QE for selected t p h 
and experimental data from Fig. [7J (b) plotted on the unitless 
pump energy scale (i.e., normalized per corresponding band 
gap energy, E g ). Dashed lines are extrapolation of the solid 
lines linear regions to determine the AET values (see Table [IJ. 



teraction is determined by the volume prefactor (v/V) 2 
which does not enter the II rate, and therefore, has no 
effect on the rate. 

In contrast, the volume prefactor appears in the ex- 
pression for the Auger recombination rate (Eq. (22 1), and 
as we show in the inset to Fig. [6j makes this effect negli- 
gible. On the frequency scale, the Auger recombination 



rate drops as Icar 



further decreasing its contribu- 



tion to QE at high energies. As a result, the only region 
where the Auger recombination processes are significant 
is near the band edge. However, this region has no sig- 
nificant contribution to QE, and we conclude that the 
Auger recombination processes can be fully neglected in 
the present analysis. 

Since the Auger recombination processes are weak, ef- 
ficient CM during the population relaxation is expected if 
the II rates are comparable to or higher than the phonon- 
assisted population relaxation rates. To evaluate the up- 
per boundary for QE in NCs, we first calculate total 
QE in bulk PbSe for typical population relaxation times 
Tph = 0.5, 1.0, and 2 ps. The results are plotted on the 
pump energy scale in Fig. [jj (a) and also compared with 
published experimental data. In general, the calculations 
reproduce the experimental trends and provide the best 
fit at Tph = 1.0 ps. The discrepancy between theory and 
experiment in the interval 1.0 < huj prn < 3.5 eV is due to 
the phenomenological nature of the phonon-assisted re- 
laxation model that lacks exact knowledge on the spectral 
dependence of the electron-phonon coupling and phonon 
DOS. Above 3.5 eV the effect of triexciton generation and 
possibly the contribution from the higher energy bands 
take place. Hence, our theory is valid below this energy 



value. 

For comparison, the contribution of the photogener- 
ation to the total QE is shown in Fig. [jj (a). Its small 
contribution is explained by the fact that the photogener- 
ation pathways reduce to a single II event taking place on 
the short (subpicosecond) dephasing timescale (Eq. (30)). 
In contrast, the population relaxation dynamics occurs 
on a much longer (picosecond) timescale allowing for the 
multiple II events that make major input into the total 
QE. 

Comparison of the calculated total QE in both NC 
and bulk PbSe with experimental result#2E2 is shown 
in Fig. [jj (b). The total QE dependence on the NCs 
diameter both in theory and experiment follows the same 
trend as the II rate. Solid lines in Fig. [jj (b), show QE 
calculated for T p h = 1.0 ps which gives the best fit for 
the bulk. However, taking into account that the upper 
values of the measured QE in the NCs are associated 
with the NCs diameter of d = 8 nm, we had to increase 
the phonon relaxation time up to r p h — 2 ps for the 
NC with d — 8 nm (blue dash) in order to get better 
agreement with experiment. The observed increase of T p h 
can be rationalized by the quantum-confinement-induced 
increase of the level spacing.^ 

Determination of AET is more convenient to perform 
using the unitless energy scale, i.e., the pump energy nor- 
malized per corresponding E g of NC or bulk. Fig. M 
presents the same curves as in Fig. [jj (b) for selected 
T p h plotted on this scale. Dashed lines extrapolate the 
linear portion of the curves at the energies below 3.5 eV. 
Their intercept points with the horizontal, QE = 1, line 
provide the AET values. The deviation from linear be- 
havior at low energy ends and apparent CM below the en- 
ergy conservation threshold, 2E g (red and green curves) 
are merely due to the ensemble averaging entering our 
calculations!^ 

The calculated values of AET are summarized in Ta- 
ble [TJ For the NCs with d = 5 nm and d = 8 nm, 
the calculated AETw 2.2E g is in excellent agreement 
with experiment.^ Furthermore, the photogeneration 
AE T=1.5 eV for d = 15 nm PbSe NCs (determined in 
Sec . |IHA| Fig. [I} is close to that, 1.4 eV, listed in Ta- 
blc[l| The latt er is an additional indication of the "pure" 
II nature of CM. For the bulk however, there is a disagree- 
ment between the calculated, 3. 7E g , and experimentally 
observed, GE g , AETpfil The discrepancy is most likely 
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TABLE I: AET in PbSe NCs and bulk calculated as the in- 
tercept point between dashed lines and the solid horizontal 
line at QE=1 in Fig. [8] 
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FIG. 9: Calculated dependence of the total QE for PbSe NC 
of d = 8 nm on the pump pulse duration r pm . The relaxation 
time is set to t p h — 1 ps. 



related to the phenomenological model of the phonon- 
assisted relaxation. Getting back to Fig. [7] (a) , one can 
clearly see that the best fit to the initial rise of QE can 
be obtained for T p h = 0.5 ps. However, for adopted 
Tph = 1 ps the model fails to accurately reproduce the ini- 
tial slope resulting in the underestimated values of AET. 

All calculations discussed above are performed for the 
cw pulses typically used in ultrafast spectroscopic stud- 
ies. However, estimated solar light correlation time is 
about 5 fs. Hence in light of photovoltaic applications, it 
is natural to calculate QE as a function of the pump pulse 
duration, r pm . Such a dependence evaluated for various 
photon energies, Hoj pm , is shown in Fig. [9] First of all, 
we notice that the cw regime is reached at r pm > 10 fs. 
For r pm < 10 fs, the QE associated with frw pm = 3.15 eV 
shows no variation whereas QE associated with the lower 
excitation energies show insignificantly small increase. 
Observed weak dependence of QE on the pulse duration 
suggests that, for the gaussian pulses, the expected in- 
crease in the QE due to the increase in the pulse duration 
is equally compensated by the reduction in the number 
of states populated by the pulse whose spectral widths is 
narrowed. 



IV. COMPARISON OF QE IN PbS AND PbSe 

Both PbS and PbSe have electronic structures de- 
scribed by the same KW effective mass model with some 
different parameters.^ This implies that our conclusions 
on the interplay between different CM pathways in PbSe 
fully apply to PbS NCs and bulk materials. In this sec- 
tion, we perform comparison of the key quantities such as 
exciton/biexciton DOS, effective Coulomb terms, the II 
rates, and the resulting QEs calculated for these semicon- 
ductors. Our numerical tests show that these quantities 
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FIG. 10: Energy dependence of (a) exciton DOS, (b) biex- 
citon DOS, and (c) effective Coulomb term in the bulk PbS 
(solid black) and PbSe (dashed black). In all panels, red lines 
show corresponding quantities calculated using PbSe parame- 
ters except for the band gap energy replaced with that of PbS. 
In panel (c), green line shows the effective Coulomb term cal- 
culated for PbSe with both band gap energy and the effective 
masses substituted with their corresponding PbS values. 



are weakly affected by the differences in the Kane mo- 
mentum and the unit cell volume. Hence, the dominant 
contributions come from the interplay between the band 
gap energies, carriers effective masses, and dielectric con- 
stants. 



Panels (a) and (b) in Fig. 10 clearly show that both 
the exciton (Eqs. (Ill])) and biexciton (Eq. Q) DOS 
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FIG. 11: Energy dependence of the photogeneration QE in the 
bulk PbS (solid black) and PbSe (dashed black) . Red (green) 
line shows the QE calculated for PbSe with the band gap 
energy (band gap energy and the effective masses) replaced 
with that (those) of PbS. 



in bulk PbS exceed the corresponding DOS in PbSe 
across the whole spectral range of interest. Compared 
to PbS {Eg = 0.28 eV), PbSe has larger band gap energy 
{Eg = 0.41 eV) and heavier carriers effective masses.^ 
In order to isolate the band gap effect, we plot the PbSe 
DOS calculated with the band gap energy replaced by 
that of PbS (red dash). This replacement results in a 
reduction of the DOS values (red dash). Hence, the ma- 
jor contribution to the steeper growth of the exciton and 
biexciton DOS in PbS comes from the heavier effective 
masses. 
Fig. 



V. 



10 (c) compares the effective Coulomb terms, 



eff 



(Eq. (151) calculated for bulk PbSe and PbS. In 
contrast to tne DOS, the effective Coulomb interaction 
in PbS is weaker than that in PbSe. To understand the 
trends, we alter the PbSe band gap value in the same way 

rx,xa 
eff 



as for the DOS, and observe small reduction in V 1 ' 1 '' 



(red dash). Altering the effective masses of PbSe with 
those of PbS, significantly reduces the V®ff X lowering 
the curve (green dash) below V**** in PbS (solid black). 
The gap between these curves is purely due to the dif- 
ference in the dielectric constants which in PbS (e = 17) 
is lower than in PbSe (e = 23). Therefore according to 
Fig. 10 (c), the narrower band gap value and lighter ef- 



rX,xx 

'eff 



making it smaller than in PbSe. 

Next, we compare the photogeneration QE in bulk PbS 
and PbSe plotted in Fig. [IT] The plot also shows photo- 
generation QE values for PbSe with the band gap energy 
and effective masses altered in the same way as in the case 
of the effective Coulomb interaction. (The color code is 
the same as in Fig. 



10 



(c).) By taking into account the 
linear dependence of the photogeneration QE on the II 
rate (Eq. (30)), one can conclude that the trends in the 
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FIG. 12: Energy dependence of (a) II rate and (b) total QE 
calculated for PbS (solid lines) and PbSe (dash) NCs, and the 
bulk. For all curves in panel (b), the phonon relaxation time 
is set to r p h = 1.0 ps. 



behavior of the QE is a result of the interplay between 
the trends in the biexciton DOS and effective Coulomb 
term shown in panels (b) and (c) of Fig. 10 respec- 
tively. Specifically, the lower dielectric constant of PbS 
contributes to the increase in the photogeneration QE 
and the competing effect of the lighter effective masses 
on the DOS and results in additional increase of 

the QE in bulk PbS. This net increase fully overruns the 
reductions associated with the band gap decrease making 
the photogeneration QE and the II rate (Fig. 
bulk PbS to be higher than in bulk PbSe. 

To demonstrate that the size quantization does not 
change the trends, we plot the II rate calculated for PbS 
and PbSe NCs and the bulk in Fig. 12 (a). It is clear from 
the plot that the PbS II rate always exceeds the PbSe II 
for identical in diameter NCs and for the bulk. Similar 
to PbSe, the PbS II rate monotonically increases with 
the NC diameter but does not exceed the bulk values. 
Finally in panel (b) of Fig. 12 we plot the calculated 
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14 



total QE of PbS and PbSe NCs and the bulk using the 
same phonon relaxation time set to T p h = 1.0 ps. The 
trends are the same as in panel (a) demonstrating that 
the overall calculated QE in PbS is higher than in PbSe. 
The latter is a result of the interplay between the band 
gap energies, effective masses and the dielectric constant 
values as discussed above. 

Ultrafast measurements of QE in bulk PbS and PbSe, 
show that QE for both materials are approximately the 
samePS However, the atomistic calculations reported 
along with the experimental data demonstrate that for 
identical phonon relaxation times (specifically, r v h = 
0.5 ps) the calculated QE in PbS exceeds that in PbSe. 
This observation is in direct agreement with our calcu- 
lations. Taking into account that both models show the 
same trend in QE for identical phonon relaxation times, 
we conclude that the II rate in bulk PbS should be higher 
than in PbSe, and consequently, the photogeneration QE. 
Furthermore, we extend this conclusion to PbS NCs ac- 
cording to the results shown in Fig. [12] 

Recently reported ultrafast measurements demon- 
strated significant reduction in QE of large PbS NCs com- 
pared to PbSe NCsP^The authors rationalize this obser- 
vation by estimating the phonon-induced energy loss rate 
which in PbS turns out to be factor two faster compared 
to PbSe. Rigorous calculation of the phonon-induced re- 
laxation is beyond our model capability forcing us to 
introduce phenomenological r p h- Since remaining elec- 
tronic structure parameters in our and reported atom- 
istic models are well validated, the calculations directly 
support the idea of the QE reduction in PbS merely due 
to the fast phonon-induced relaxation processes charac- 
terized by Tph- 



V. DISCUSSION AND CONCLUSIONS 

The effective mass model adopted in this paper has 
significant limitations accounting only for the L-valley 
optical transitions.^ It fails to catch contributions from 
the higher energy transitions originating at other Bril- 
louin zone points. Specifically in PbSe, the S-point tran- 
sitions show large contribution to the absorption spec- 
tra above 1.6 eV.^ Experimentally, the CM dynamics 
is probed using ultrafast transient absorption and fluo- 
rescence techniques.^ These techniques explicitly deter- 
mine the number of carriers accumulated near the band 
edge of the lowest in energy L-valley. Following the ex- 
citation with a high energy optical pump, the excitons 
and biexcitons can be generated not only in L-valley (de- 
noted L-excitons and biexcitons) but also at S-point (de- 
noted S-excitons and biexcitons). The contribution of 
the E-excitons to the transient ultrafast signal depends 
on the mixing between the L and S-points. Rigorous 
calculations of the mixing can be done through atom- 
istic calculations. However, such calculations become 
tremendously expensive for the considered large diameter 
NCs. Below, we provide quantitative discussion of possi- 



ble mixing mechanisms and rough estimates of the their 
contribution to QE for the pump energy below 3.5 eV. 

First, we consider mixing between valleys due to the 
size quantization which could be efficient in very small 
(d < 3 nm) NCs at low excitation energies. How- 
ever, for the large NCs, considered here, the excitations 
with the ene rgies higher than 2.5 eV are effectively of 
the bulk-typ e! 76 * 77 ^ resulting in the negligible quantum- 
confinement-induced mixing. Another option is the II 
event initiated by the optically prepared E-exciton and 
resulting in the direct production of an extra L-exciton. 
Specifically, the excess energy of the S-exciton is trans- 
ferred to create an extra electron-hole pair within L- 
valley. The AET for this process can be roughly esti- 
mated as AET^ L « AET L + (Ef - Eg) = 3.3 eV point- 
ing to an efficient CM at photon energies hu) pm > 4 eV. 
Note that efficient CM events fully constrained to E-point 
should occur at even higher energies which are above the 
photon energy range considered in this paper. 

Finally, the optically prepared S-exciton can partic- 
ipate in the II event resulting in the creation of the 
electron-hole pair through the inter- valley transition. For 
instance, a higher energy conduction band electron from 
S-point can be transferred to L-valley releasing in the 
excess energy to create an additional electron-hole pair 
through the L-valence band to E-conduction band tran- 
sition. As we mentioned above, most of the states in- 
volved in efficient CM are of the bulk-type with the quasi- 
momentum being a "good" quantum number. Therefore, 
the Coulomb matrix elements describing the II processes 
become ~ 1/| Aq| 2 where Aq is the difference between the 
electron/hole initial and final states. Accordingly, the II 
events involving the inter-valley transitions require large 
variation in the quasi-momcntum, Aq, and, therefore, 
becomes less favorable. 

Experimental verification of the discussed interplay be- 
tween the CM pathways contributing to the photogen- 
eration QE requires direct measurements of the inter- 
band Coulomb interactions. Distinguishing the inter- 
and intraband Coulomb interactions is a challenging task, 
since in the transient absorption and time-resolved fluo- 
rescence experiments these two components contribute 
to a measured total energy shift between the exciton 
and biexciton bands. Recently, we have theoretically 
shown that the two-dimensional double-quantum coher- 
ence spectroscopy is capable to probe directly the inter- 
band Coulomb interactions.^ Alternatively, the photo- 
generation QE can be determined based on its linear re- 
lationship to the II rate which has been established above 
(Eq. ([30])). The later rate can be obtained by fitting the 
total QE as a function of the excitation energy™ How- 
ever, the determination of the effective dephasing rate, 
also entering Eq. ( 30 ) , might require use of coherent ex- 



perimental technique (e.g., photon echo) and/or addi- 
tional theoretical calculations. 

The absolute photon energy scale, hw pm , is used 
throughout this paper (except Fig. ^ a) and [8). This 
scale is more useful to discuss fundamental physical 
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mechanisms behind t he CM dynamics, e.g., the role of the 
quantum confinement 1251231 However, performance of pho- 
tovoltaic devices can be ch aracterized by various power 
conversion efficiencies! 34 * 63 ! In particular, by the QE cal- 
culated on a dimensionless energy scale which is nor- 
malized per NC or bulk band gap, hLu pm /E g . Despite 
lower values of QE in NCs compared to the bulk which 
show up on the absolute photon energy scale, the perfor- 
mance of prospective photovoltaic d evices b ased on NCs 
can overrun their bulk counterpart! 35 * 63 ^ By compar- 
ing Fig. H] with Fig. [7| (b), one can notice that smaller 
size NCs have better rise in QE than large ones and 
the bulk. This results from the quantum-confinement- 
induced blue-shift of the band gap energy.^ In light of 
this effect and the requirement to match solar radiation 
peak energy, semimetal NCs can become more efficient 
solar energy converters.^ 

Further enhancement of the interband Coulomb inter- 
actions in NCs can possibly be reached in NC-metal het- 
erostructures in which broad surface plasmon response 
is tuned in resonance with the exciton states participat- 
ing in CM. This issue requires additional study, since 
besides the enhancement of the Coulomb interactions, 
the surface-plasmons- induced ohmic and radiative energy 
losses can become significant!^ On the other hand, ad- 
ditional enhancement of QE can in principle be achieved 
by adding the surface states and ligands. This question is 
difficult to analyze using adopted effective mass appro x- 
imation. However, the atomistic calculation s! 30 * 31 * 62 * 82 ! ^ 
combined with the IESM should be extremely helpful in 
addressing this problem. 

In conclusion, we have performed systematic numeri- 
cal investigation of CM mechanisms in NC and bulk PbSe 
and PbS using our IESM parametrized by the effective 
mass KW model. We focused on the role of quantum 
confinement on the photogeneration and total QE. The 
analysis of the photogeneration pathways resulted in un- 
expected conclusion that the photogeneration processes 
reduce to a single II event due to the complete destruction 
of the associated pathways interference. This allowed us 
to explains the minor role of the photogeneration in total 
QE dominated by the multiple II events occurring during 
the phonon-induced population cooling. Comparison of 
the size scaling of the effective Coulomb interaction and 
the biexciton DOS in transition from NC to the bulk limit 
showed that weak enhancement of the former quantity is 
overrun by a significant reduction in the latter one. This 
explains the higher values of the QE in bulk compared 
to NCs as plotted on the absolute energy scale, and well 
agrees with the previously reported experimental studies 
and some theoretical predictions. However, the quantum 
confinement induced increase in E g makes NCs more ef- 
ficient than bulk for practical photovoltaic applications. 
We have also found that the variation in the pump pulse 
duration does not significantly change the QE. Compar- 
ison of QE in PbSc and PbS suggests that the II pro- 
cesses are more efficient in PbS. However, variation of 
the material-dependent cooling time can strongly affect 



the total QE. Finally, we have identified the limitations 
of our model and defined its applicability range. 
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Appendix A: Volume scalings of DOS, transition 
dipoles, and the effective Coulomb interaction 

In this Appendix, we derive the volume normalization 
prefactors entering Eqs. (11 1— ( 16 ) . The prefactors can- 



cel out the volume dependence of the latter quantities 
making them intensive variables in the bulk (i.e., ther- 
modynamic) limit. To preserve the dimensionality of the 
intensive variables, we use the V/v ratio instead of V 
where v is the unit cell volume. This ratio defines number 
of unit cells and goes to infinity in the bulk limit. Here, 
we also use the relationships derived in Appendix [5] that 
connect the transition dipole movements and Coulomb 
matrix elements represented in the single-particle (KW) 
carrier basis set and in the exciton/biexciton basis. 

We start with the simple fact that in the bulk limit, the 
single-particle DOS is proportional to the system volume 
VP^Defining the exciton (biexciton) DOS as a joint DOS 
of an electron and a hole (two electrons and two holes), 
one immediately finds that they have V 2 (V 4 ) scalings. 
Therefore, we introduced v 2 /V 2 (v A /V A ) prefactor into 
Eq. {11) (Eq. Q). 

To obtain the volume prefactors in the optically al- 



lowed DOS given by Eqs. (13) and (14 1, we first notice 



that the transition dipole matrix element for a single car- 
rier has no volume scaling and accounts for the total mo- 
mentum conservation, i.e., Mij ~ V^M^^S^ ^ .^Fur- 



ther using Eq. (B2), we find that the optically allowed 



exciton DOS scales as 



X>Sol 2 % (ai) 

a 



Here and below, we use the same argument as in the DOS 
analysis that for a fixed energy interval, J2 kp ~ V. Using 
the first term in the expansion of the intraband biexciton 



dipole moment given by Eq. (B3 ), one finds that the joint 
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optically allowed biexciton DOS scales as 

E |/4fT^i - u%°)8(w2 - wf) ~ (A2) 

E E i^Ckj^^i - <4„ - < - < - O 



x<5(w 2 - wj< 



According to Eqs. ( Al ) and ( A2 1, we introduced the pref- 
actors v/V and (v/V) 4 into Eqs. (13) and (14), respec- 
tively. 

To determine the volume scaling of the effective 
Coulomb term (Eq. (15)), we first evaluate the scaling 



of the following auxiliary quantity 
£( Wl> lo 2 ) = E IKX I 2 <^i - <) 6 ^ - O- ( A3 ) 



According to Eq. (BIO), the Coulomb matrix ele- 
ments in the free carriers basis scales as Vj^ki ~ 



V 1 VSt i k i ,k*k,5k 4 -k l) k J -k Jk - Then using Eq. <|B6j), one 
finds that 



K(-^2)~r 2 E l^, kfc , kl+kj+k J 2 (A4) 

kik J kfek i 

x<5(wi-a; k . +k ._ kfe -^)-l/ 2 . 

By taking into account that /C should be normalized by 
the exciton and biexciton DOS and tak e sq uare root, we 
introduce the prefactor (V/v) 2 into Eq. ( 15 ). The volume 
prefactor in Eq. ( 16 ) can be obtained in the same way 



with the help of Eqs. (B5|, (BIO), and (BID 



Appendix B: Representation of transition dipoles 
and interband Coulomb matrix elements in KW 
basis set 

In this Appendix, we provide closed expressions for 
the exciton and biexciton transition dipoles and the in- 
terband Coulomb matrix elements used in the numerical 
calculations. To derive these expressions we follow the 
procedure outlined in Appendices A-C of Ref. [57]. 

The second quantization is performed using the basis 
of KW states, {^(r)} defined in Eq. (26), by introducing 
the following field operators 

*(r) = E [e(Ei)*i(T)ci + B(-Ei)*i( r )4] '( B1 ) 

i 
i 

where Q(E) is the step function, Cj and di (cj and d\) 
are conduction band electron and valence band hole an- 
nihilation (creation) operators, respectively. 



Using this representation and the definition of the ex- 
citon and biexciton states given by Eqs. (27) and (28), 



it is straight forward to show that the transition matrix 
elements entering Eqs. ^ and Q are 



XX 



(x a |M|x ) = 
[xx k \M\xxi] 



[Mp p 5 qq > 
[Sppf Sqqf — 

[M h JlS rr , 



±y± qr > 

M e q ? p 6 qp , - 



5 rs '8 r > 

MpgSpq, 



(B2) 
(B3) 



Spq' $p'q\ 
rhh , 



M n / r 5 sr , - M^ s 5 rs , 



+ M q>q$pp' 



respectively. In Eq. (B2), the exciton index a = {<?, r} 



and is the interband electron-hole transition dipolc 
matrix element calculated in the KW basis set. In 
Eq. (B3), the biexciton indices are k — {p'q',r's'} and 
I = {pq, rs} as well as MfJ and M^ 1 are the matrix 
elements of the intraband transition dipole operator cal- 
culated using the electron and hole KW wave functions, 
respectively. 

In the envelope function approximation, the KW tran- 
sition matrix elements entering Eqs. (B2) and (B3) have 
the following generic forrrPl 



(2) 



Mg> = E 

m— 1 



dv KWl'p^fr), 
Mg> = P t z J dv (\Fl{v)]*Fi(v) + [^ ; (r)]*F/(r) 



-[Fi(v)]*Fl(r)-[Fi(r) 



*FZ( 



(B4) 



where the indices i,j can be associated with the elec- 
tron and hole states, and p = — iSV. Pi stands for the 
longitudinal dipolc moment component of the band-edge 
Bloch function, u m (r), and z is the unit vector in (111) 
direction of the PbSe (PbS) lattice. Both M« and M^ 2 ) 
can be evaluated analytically in the bulk limit, where 
M^ 1 ' vanishes identically.^ 

Using the same approach as above, one can show that 
the Coulo mb matrix elements entering Eqs. ([5|, ( fl5| , 
(fl6|, and (|A4l are 



■xrxxfi 
V l,0 



k . a 



(x Xl \v\x ) = v;^ - v; q f r \ (B5) 

( XXk \V\x a ) = (V q %^q - Vffiq) 6 s ,, r 

i { -ireeeh yeee/t \ r 

-|- yvq'p'qs' Vp'q'qs 1 ) u r' ,r 

I ( -rrehhh _ vrehhh \ c 
T" \ V q' rs ' r > Vq'rr's') ^p' ,Q 

I / jrehhh -irehhh \ r 

' v p'rr's' p'rr's') Q f ,Q m 

These quantities depend on the long-range contributions 
whose matrix elements in the KW basis read 



(B6) 



Vijki — 



dvidr 2 



qri - r 2 



4 

E 

rn — 1 



1 rn—1 

xV[i*(r a )]*F*(r 2 ). (B7) 
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Here, e denotes the screened dielectric function values 
evaluated at the optical frequencies. For bulk PbSe and 
PbS, we set e bulk =23 and e bulk =17, respectively.^ The 
dielectric constant in the NCs has been evaluated using 
the following expression^ 



JVC 



(d) = 1 + (6 



bulk 



(£ bulk 



AE) 2 



AE) 2 ' 



(B8) 



where E^ nlk + AE = 2.73 eV and 3.14 eV are the energy 
of the first pronounce d abso rption peak in the bulk PbSc 
and PbS, respectively! 221221 

To derive the volume scaling of the Coulomb matrix 
elements given by Eq. (B7), we assume the bulk limit 



in which the envelope functions become plane waves, i.e. 
F4(ri) = F^e^'" 1 /VV. In this basis, the Coulomb ma- 
trix elements can be written as 



V--ui - — — VlF F 3 ] * F l F k 

1 



I I dr 1 dr 2] 



(B9) 

ikiri— ikjk2+ik fc r 2 +ik i ri 



The integral evaluation leads to the final expression 

Vijkl = yz ^k i ,k :f ,k )fe ,k t ^k i -k,,k ;j -k fc , (BIO) 



where 



ki,k,,k fe ,ki 



e2 V- 



\pi pj]* pi pk_ 
L m n\ m n i 



4tt 



(Bll) 



Eqs. (BIO) and (Bll) clearly demonstrate a general prop- 



erty that the Coulomb matrix elements used in the nu- 
merical calculations of the bulk limit scale inversely pro- 
portional to the volumePS 

The nonparabolicity in the employed k • p-Hamiltonian 
is crucial for the accurate evaluation of the above 
Coulomb matrix elements. Specifically, summations over 
the spinor components in Eq. (B7) implies that the in- 



terband Coulomb scattering amplitudes vanish exactly 
if the non-diagonal terms of the k • p-Hamiltonian are 
set to zero, i.e., in strictly parabolic caseP^In the bulk, 
where quasimomentum A: is a "good" quantum number, 
diagonal and off-diagonal matrix elements of the k • p- 
Hamiltonian scale as k 2 and k, respectively. At high 
energies, where the diagonal elements dominate over the 
off-diagonal ones, the latter can be treated perturbatively 



giving k 



-1/2 



as the contribution of the hole (elec- 



tron) states to a high energy electron (hole) wave func- 
tion. 
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